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We consider the effect of rescattering of pairs of quasiparticles in the Cooper channel resulting in 
the strong renormalization of second-order corrections to the spin susceptibility in a two-dimensional 
electron system. We use the Fourier expansion of the scattering potential in the vicinity of the 
Fermi surface to find that each harmonic becomes renormalized independently. Since some of those 
harmonics are negative, the first derivative of the spin susceptibility is bound to be negative at 
small momenta, in contrast to the lowest order perturbation theory result, which predicts a positive 
slope. We present in detail an effective method to calculate diagrammatically corrections to the 
spin susceptibility to infinite order. 
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I. INTRODUCTION 

The study of the thermodynamic as well as micro- 
scopic properties of Fermi-liquid systems has a long 
history^i^i^i^ but the interest in nonanalytic corrections 
to the Fermi-liquid behavior is more recent. The exis- 
tence of well-defined quasiparticles at the Fermi surface 
is the basis for the phenomenological description due to 
Landau^ and justifies the fact that a system of interact- 
ing fermions is similar in many ways to the Fermi gas. 
The Landau theory of the Fermi liquid is a fundamental 
paradigm which has been successful in describing prop- 
erties of ^He, metals, and two-dimensional electronic sys- 
tems. In particular, the leading temperature dependence 
of the specific heat or the spin susceptibility (i.e., Cs 
linear in T and Xs approaching a constant) is found to 
be valid experimentally and in microscopic calculations. 
However, deviations from the ideal Fermi gas behavior 
exist in the subleading terms. 

For example, while the low-temperature dependence 
of Cs/T for a Fermi gas is a regular expansion in T^, a 
correction to Cs/T of the form InT was found in three 
dimensions.^ These nonanalytic features are enhanced in 
two dimensions and, in fact, a correction linear in T is 
foundj^i^i^ These effects were observed in ^He, both in 
the three-^ and two-dimensional caseJ^ 

The nonanalytic corrections manifest themselves not 
only in the temperature dependence. For the special case 
of the spin susceptibility, it is of particular interest to de- 
termine also its dependence on the wave vector Q. The 
deviation Sxs from the T = Q = value parallels the 
temperature dependence of the specific heat discussed 
above: from a second-order calculation in the electron in- 
teraction, corrections proportional to In Q and Q were 
obtained in three and two dimensions respectively^i^i^ 
On the other hand, the dependence on T was found to be 
6xs ^ in three dimensions^iii (without any logarith- 
mic factor) and 5xs ~ T in two dimensions. 



We cite here the final results in the two dimensional case 
(on which we focus in this paper), valid to second order 
in the interaction potential V{q), 
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Here m is the effective mass, /cf is the Fermi wave vec- 
tor, Ef = k'^/2m^ and we use h = 1 throughout the 
paper. Our purpose is to extend this perturbative re- 
sult to higher order by taking into account the Cooper 
channel renormalization of the scattering amplitudes. 

The extension to higher order of the second-order 
results has mostly focused on the temperature depen- 
dence, both for the specific heat^^^^^^^^^^^^^^ and the spin 
susceptibilityJ^i^i2^i2^ Recently the spin susceptibility 
has been measured in a silicon inversion layer as a func- 
tion of temperature.^^ A strong dependence on T is ob- 
served, seemingly incompatible with a Fermi-liquid 
correction, and the measurements also reveal that the 
(positive) value of the spin susceptibility is decreasing 
with temperature, in disagreement with the lowest or- 
der result cited above. This discrepancy has stimulated 
further theoretical investigations in the nonperturbative 
regime. Possible mechanisms that lead to a negative 
slope were proposed in Refs. 21 and 22 if strong renor- 
malization effects in the Cooper channel become impor- 
tant. These can drastically change the picture given by 
the lowest order perturbation theory, allowing for a non- 
monotonic behavior and, in particular, a negative slope 
at small temperatures. 
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The mechanism we consider here to modify the hnear 
Q dependence is very much related to Ref. [2l|. There it 
is found that, at Q = and finite temperature, y^(2/ci?) 
in Eq. ([T]) is substituted by |r(7r)p, where 



Pa* + Qp. 



(4) 



is the scattering amphtude in the Cooper channel with 
9 being the scattering angle = tt corresponds to the 
backscattering process). An additional temperature de- 
pendence arises from the renormalization of the Fourier 
amplitudes 



Tn{kBT) 



Vn 
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where is a large energy scale W ^ Ep and Vn are 
the Fourier amplitudes of the interaction potential for 
scattering in the vicinity of the Fermi surface 



V{2kF sme/2) = 
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A negative slope of Sxs is possible, for sufficiently small 
T if one of the amplitudes Vn is negative j^^i^^ ^^ For 
^In^^^ = 1, the denominator in Eq. © di- 
verges what corresponds to the Kohn-Luttinger (KL) 
instability.^^ At T > Tkl the derivative of the spin sus- 
ceptibility is negative due to the singularity in VniksT) 
and becomes positive far away from Tkl- 

At T = an analogous effect occurs for the momen- 
tum dependence. Indeed, it is widely expected that the 
functional form of the spin susceptibility in terms of ksT 
or vpQ is similar. As in the case of a finite temperature, 
the lowest order expression gains an additional nontriv- 
ial dependence on Q due to the renormalization of the 
backscattering amplitude y^(2/ci?). We obtain 



6xs{Q)=2K{Q,Q)\T{Trf 
where r(7r) is given by Eq. @ and 
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Such result is obtained from renormalization of the in- 
teraction in the Cooper channel, while other possible ef- 
fects are neglected. Moreover, at each perturbative or- 
der, only the leading term in the limit of small Q is kept. 
Therefore, corrections to Eq. ([7j) exist which, for exam- 
ple, would modify the proportionality of 5xs to |r(7r)p 
(see Ref. 21). However, in the region vpQ ^ ^bTkl^ 
close to the divergence of TnivpQ) relative to the most 
negative Kz, Eq. ([7j) is expected to give the most impor- 
tant contribution to the spin susceptibility. 

The result of Eqs. ([7j) and ([8|) could have been per- 
haps easily anticipated and, in fact, it was suggested 
already in Ref. [23. The question of the functional de- 
pendence of the spin susceptibility on momentum is cru- 
cial in light of the ongoing studies on the nuclear spin 



FIG. 1: The building block (on the left) of any ladder diagram 
(on the right). Of special interest is the limit of correlated 
momenta p = — k, leading to the Cooper instability. 



ferromagnetism^2L2^i2S as the stability of the ferromag- 
netic phase is governed by the electron spin susceptibil- 
ity. In this context, Eqs. ([7j) and ([8|) were motivated 
by a renormalization-group argument. We provide here 
a complete derivation, based on the standard diagram- 
matic approach. 

The paper is organized as follows: in Sec. [Ill we discuss 
the origin of Cooper instability and derive expressions for 
a general ladder diagram, which is an essential ingredient 
for the higher order corrections to the spin susceptibility. 
In Sec. mil we give a short overview of the lowest order 
results to understand the origin of the nonanalytic cor- 
rections. Based on the results of Section [III we provide an 
alternative derivation of one of the contributions, which 
can be easily generalized to higher order. Sec. II VI contains 
the main finding of this paper: the Cooper renormaliza- 
tion of the nonanalytic correction to the spin suscepti- 
bility is obtained there. We find an efficient approach 
to calculate higher order diagrams based on the second- 
order result. In Section [Vj the diagrammatic calculation 
is discussed in relation to the renormalization-group ar- 
gument of Ref. 27. Sec. I VII contains our concluding re- 
marks. More technical details have been moved to the 
Appendixes [Al[Dl 



II. PARTICLE-PARTICLE PROPAGATOR 

In this section we consider a generic particle-particle 
propagator, which includes n interaction lines, as de- 
picted in Fig. [H The incoming and outgoing frequen- 



cies and momenta are k 



and k'^.p'^, 



respectively, us- 
ing the relativistic notation k^ — (cc;/e,k). This particle- 
particle propagator represents an essential part of the di- 
agrams considered in this paper and corresponds to the 
following expression: 

n-l 

X n ^(^^ - + ft,M)^(|qi+i - qil), (9) 

where = — p. The frequencies are along the imag- 
inary axis, i.e., G{k^) = G(a;/c,k) = {icok — ^k)""*^, where 
k'^/2m - Ef with A: = |k|. 
In particular, we are interested in the case when the 
sum of incoming frequencies and momenta is small; i.e.. 



0. Under this assumption we obtain the 



following useful result for which we provide details of the 
derivation in Appendix {K\ 

n(") {P^,9)= Y.' "m, (Pp) . . . nM„_i (P^) 

Mi...M„_i 

X V^M,...M„_,(^P,^), (10) 

where the sum is restricted to Mi = 0, ±2, ±4 . . .. The 
angle of P = p + k is from the direction of the incoming 
momentum p, i.e., Op = Z(P,p), while = Z(p',p). In 
the above formula, 



no(P^) = £ln 



w 



and (M even) 
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with W ~ Ef a high energy cutoff and (j) = arctan j^;^. 
Notice that IIm{Pij,) has no angular (6>p, 6) dependence, 
which is only determined by the following quantity: 



X e 



m-Mi-...Mn 



im Up—zmU r /i q\ 

0Mi^M2 + ...Mrr,-i,m' U'jj 



defined in terms of the amplitudes Vn- Equation (|6]) 
can be used to approximate the interaction potential in 
Eq. (|9]) since the relevant contribution originates from the 
region of external (p ^ p' ^ k ^ k' ^ hp) and internal 
momenta (|p + qi| ^ |k — q^| ^ kp) close to the Fermi 
surface. Furthermore, the direction of P can be equiva- 
lently measured from k without affecting the result since 
Op = Z(P, k) + TT and e^^'^ = 1 {m' is even). 



Notice also that the leading contribution to Eq. ([TO]), 
in the limit of small Qp and P, is determined by the 
standard logarithmic singularity of no(P^). However, 
it will become apparent that this leading contribution 
is not sufficient to obtain the correct result for the de- 
sired (linear-in-Q) corrections to the response function. 
The remaining terms, IIm{Pi^), are important because 
of their nonanalytic form due to the dependence on the 



ratio 



KM 
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III. 



SECOND-ORDER CALCULATION 



The lowest order nonanalytic correction to the spin 
susceptibility has been calculated in Ref. ^ as a sum of 
four distinct contributions from the diagrams in Fig. [21 



V\0)], (14) 
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6x^^\Q) = K{0,Q)V{0)V{2kp), (16) 



Sx^^\Q)=K{0,Q)[V\2kp) 




FIG. 2: The nonvanishing second-order diagrams contributing 
to the nonanalytic behavior of the electron spin susceptibility. 




FIG. 3: Labeling of the 5xi^ dia gram, as in Eq. (|18p. 



(2) (2) 

and — ~^x\ such that the final result reads as 



5xf\Q)=2K{Q,Q)V\2kF). 



(17) 



We refer to Ref. 8 for a thorough discussion of these 
lowest order results, but we find it useful to reproduce 
here the result for 5x^i^ - In fact, Eq. ([T4j) has been ob- 
tained in Ref. 8 as a sum of two nonanalytic contributions 
from the particle-hole bubble at small {q = 0) and large 
{q = 2kp) momentum transfer. These two contributions, 
proportional to V^(0) and V^(2/cf), respectively, can be 
directly seen in Eq. ([T4|) . However, it is more natural for 
our purposes to obtain the same result in the particle- 
particle channel by making use of the propagator dis- 
cussed in Sec. [Ill This approach is more cumbersome 
but produces these two contributions at the same time. 
Furthermore, once the origin of the lowest order non- 
analytic correction is understood in the particle-particle 
channel, higher order results are most easily obtained. 

We start with the analytic expression of ^Xi iQ) (see 
Fig. [3]) in terms of H^^) , the n = 2 case of Eq. (fTU|) : 



r d^k r 
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It is convenient to define the angle of k as ^/^ = Z(k, Q), 
and Op = Z(P,k). We first perform the integration in 
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d^k^ as explained in Appendix \E\ to obtain 



(2) 
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A/(fJp + WpPcosOpY + {vfQY \ ,„ FIG. 4: The series of diagrams contributing to ^Xi(Q)- 

-'- 7^ '■ 7i dOp. (19) 



Following the method of Ref. 481, we rescale the integra- 
tion variables: Qp = RvpQ smc/)^ P = RQcoscj)^ and 
dQpdP = RvpQ'^dRdct). This gives 

Sx?^ = -^ R'dR #/ n(2)(i?,0,^p,o) 

TT^^F Jo Jo Jo 

. / , A/i?2(sin^ + z cos cos 6>p)2 + 1 \ 7/, 

X COS0 1 — — ■ ■ — dOp. 

\ R{sm<p -\- i cos <p cos Op) I 

(20) 

where, from Eqs. ([TQ|) and (p!3|) , 

n(2)(i?,,/,, ep,e) = ^V^(^p, 0)nM(i?, 

M 

=^'nM(J?,'/')E^™^— Me'""'""™' (21) 

M m 

with the primed sum restricted to even values of M. 

Now we can see clearly that the linear dependence on 
Q in Eq. ([2Q|) can only be modified by the presence of 
n^^^ in the integrand because of 



Uo{RA) = -\n^ 



m 
2^ 



lni?(l + sin^). (22) 



The first logarithmic term is diverging at small Q but 
does not contribute to the final result since it does not 
depend on 6>p and (j). In fact, if we keep only the 
^ In contribution, after the change of variable r = 
/?(sin^+icos0cos^p) in Eq. ([2Q|) , we obtain the angular 

integral j'^^ d O p JJ^^^ cos 0(sin cos cos Op)~^ dcj) — 
[cf. Eq. ()C6p for M = 0]. Details of the calculation are 
provided in Appendix O 

Therefore, only the second term of Eq. ([22|) is relevant. 
The integral in Eq. ([2Q|) becomes independent of Q and 
gives only a numerical prefactor. The final result is given 
by Eq. ([T4|) , in agreement with Ref. In a similar way, 
the remaining diagrams of Fig. [2] can be calculated. 



IV. HIGHER ORDER DIAGRAMS 

In this section we aim to find the renormalization of 
the four diagrams depicted in Fig. [2] due to higher order 
contributions in the particle-particle channel. It is well 
known that the scattering of two electrons with opposite 
momenta, in the presence of the Fermi sea, leads to the 
emergence of a logarithmic singularity. ^^^^^ Furthermore, 
in two dimensions there are just two processes that con- 
tribute to 6x1 \Q)^ namely, forward- (small momentum 





FIG. 5: An example of diagram contributing to 5x2{Q)- The 
maximally crossed diagram (left) is topologically equivalent to 
its untwisted counterpart (right) in which the particle-particle 
ladder appears explicitly. 



transfer, ^ = 0) and back-scattering (large momentum 
transfer, q = 2A:p). This results in the renormalization of 
the scattering amplitudes appearing in the second-order 
results (see Sec. H]). 

A direct calculation of the particle-particle propaga- 
tors, depicted in Fig. (TJ shows that for n + 1 interaction 
lines, the divergence always appears as the nth power 
of a logarithm. At each order of the perturbative ex- 
pansion, we only consider the single diagram which con- 
tributes to the nonanalytic correction with the leading 
logarithmic singularity. This requirement restricts the 
freedom of adding interaction lines in unfettered manner 
to the existing second-order diagrams: in order to pro- 
duce the most divergent logarithmic term, all interaction 
lines have to build up at most one ladder for 5xi^ ^X2, 
and 5xa^ or two ladders for Jxs- 

The subset of diagrams generated in this way is not 
sufficient to obtain the general momentum dependence 
of the spin susceptibility. However, if one of the har- 
monics Vn is negative, these diagrams are the only rele- 
vant ones in the vicinity of the Kohn-Luttinger instabil- 
ity, vpQ ^ kpTKL' Furthermore, at each order n in the 
interaction, it suffices to keep the leading contribution in 
Q of the individual diagrams. This turns out to be of or- 
der Qln"^"^ Q because the term proportional to ln"^~^ Q 
is suppressed by an additional factor Q^. Other pertur- 
bative terms, e.g., in the particle- hole channel, can be 
safely neglected as they result in logarithmic factors of 
lower order. 

In the following we discuss explicitly how to insert a 
ladder diagram into the pre-existing second-order dia- 
grams and show the line of the calculation that has to be 
carried out. 



A. Diagrams 1, 2, and 4 

These three diagrams can all be expressed to lowest or- 
der in terms of a single particle-particle propagator H^^^ , 
which at higher order is substituted by U'^^^ . For the first 
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FIG. 6: A maximally crossed diagram (left) and its untwisted 
equivalent (right) contributing to 6x4{Q)- 



term we have 



xG{k^^Q^)G{P^-k^)U^^\P^,0), (23) 

where the n = 2 case was calculated in Sec. Illli The 
corresponding diagrams are, in this case, easily identified 
and shown in Fig. [H 

(2) 

It is slightly more complicated to renormalize 

(2) 

and 6x4 . It requires one to realize that the diagrams 
depicted in Fig. [5] are topologically equivalent; i.e., the 
maximally crossed diagram on the left is equivalent to the 
untwisted ladder diagram on the right. A similar analy- 
sis shows how to lodge the ladder diagram into ^xi^^ 5 
illustrated in Fig. [6l The corresponding analytic expres- 
sions are: 



(27r)3 J (27r)3 
X G{k^ + g^)G(P^ - fc^)n(") (P^, tt), (24) 



(27r)3y (27r)3 



G(fc^)G(fc^ + Q^) 



X G{P^ - k^)G{P^ -k^- Q^)n(")(P^, tt). (25) 

We show now that the final results can be simply ob- 
tained to leading order in Q based on the second-order 
calculation. In fact, we can perform the integration in 
d^k and the rescaling of variables as before. For Sxi we 
have 



in) 



r R^dR r dop r^\^^\R,^,op,o) 

t^^vf Jo Jo Jo 



^y R^ (sm (j) -\- 1 cos (j) cos p)^ + , , , 

1 — ^^TT-. — ] ] 7r-\ cos (p d(p. 

R[sm(p -\- 1 cos (p cos Op) J 

(26) 

In the above formula, the Q dependence in the integrand 
is only due to H^^^. It is clear that a similar situation 
occurs for the second and fourth diagrams. 

The Q dependence of the rescaled Eq. (p!Q|) is deter- 
mined (as in the second order) by the factors HoiR^cj)). 
The first term appearing in no(i?, 0), see Eq. ([22|) . is 
large in the small Q limit we are interested in. There- 
fore, we can expand H^^^ in powers of ^ In and re- 
tain at each perturbative order n only the most diver- 
gent nonvanishing contribution. The detailed procedure 



is explained in Appendix [Dl It is found that the largest 
contribution from H*^^) is of order (^ In ^^)^~^. How- 
ever, as in the case of the second-order diagram discussed 
in Sec. Illlt this leading term has an analytic dependence 
on P^ (in fact, it is a constant), and gives a vanishing 
contribution to the linear-in-Q correction to the spin sus- 
ceptibility. Therefore, the (^In^^)"^"^ contribution is 
relevant here. 

A particularly useful expression is obtained upon sum- 
mation of n^^^ to infinite order. In fact, for each diagram, 
the sum of the relative series involves the particle-particle 
propagator only. Therefore, 6x1^ ^^X2, and Sxa are given 
by Eqs. ([23])-([25l) if H^^) is substituted by 



n(-)(p^,^) = ^n(-)(p^,^). 



(27) 



n=2 



The relevant contribution of 11^^^ (P^, ^), in the rescaled 
variables, is derived in Appendix [Dl The final result is 



n=2 

=Y!^m{R. <l^)Y^Tmrm-Me'^'--'^' + . . . , (28) 



M 



which should be compared directly to Eq. ([2T]) . The only 
difference is the replacement of Vn with the renormalized 
amplitudes F^, which depend on Q as in Eq. ([8|). 

Hence, it is clear that the final results follow immedi- 
ately from Eqs. ([Ti|) - (fT6l) : 



Sxi{Q)=K{0,Q)[T\0)^T\7r) 

(5x4(Q) = i^(o,Q)r(o)r(^), 



(29) 
(30) 



and (5x2 (Q) = —Sx4.{Q)' We have used notation (|4]) 
while i^(0, Q) is defined in Eq. ([3|). This explicitly proves 
what was anticipated in Sec [J (and in Ref. ^), i.e., that 
the renormalization affects only the scattering ampli- 
tude. The bare interaction potential is substituted by 
the dressed one, which incorporates the effect of other 
electrons on the scattering pair. 



B. Diagram 3 

(2) 

The last diagram Sxs differs from those already 
discussed in the sense that it allows for the separate 
renormalization of either the upper or lower interaction 
line. This results in the appearance of two equivalent 
third-order diagrams and three fourth-order diagrams (of 
which two are equal), and so forth. These lowest order 
diagrams are shown in Fig. [71 Accordingly, we define 
the quantities Sxs'^\ where ladders of order i and j are 
inserted in place of the original interaction lines. In par- 



ticular, Sx 



in) 



and 



5xm = ^xf'\Q)- 



(31) 
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FIG. 7: The series of diagrams contributing to Sx3{Q)- At 
the top, the second- and third-order order diagrams. Two 
equivalent third-order diagrams arise from the addition of a 
parallel interaction line to either the upper or the lower part of 
the second-order diagram. At the bottom, three fourth-order 
diagrams. 



The second difference stems from the fact that a finite 
nonanalytic correction is obtained from the leading terms 
in the particle-particle ladders of order (^ln^^^)*~^ 
and In ^^^)-^~^, respectively. In fact, extracting this 
leading term from Eq. (p!Q|) we obtain 



n«(P„^) = 5:i-e-»^gln^y"V... , (32) 

n 

and by performing the sum over j we get 

oo 

^n(^-)(p^,^) = r(^) + ... . (33) 

A similar argument can be repeated for the ith order in- 
teraction ladder. Therefore, the bare potential is replaced 
by renormalized expression @ and the final result, 

Sx3{Q) = K{0,Q)[T\Tr)-T\0)], (34) 

is immediately obtained from Eq. (p!5|) . 



C. Renormalized nonanalytic correction 



Combining the results of Sec . II V Al and II V Bl it is clear 
that the final result has the same form of Eq. ([T7j) if 
V{2kF) is substituted by r(7r). The explicit expression 
reads as 



SxsiQ) 



247r4 kF 



Vn{-iy 



mVn 1^ vfQ 
27r W 



(35) 



V. RELATION TO THE RG TREATMENT 

As discussed, our calculation was partially motivated 
by the renormalization group (RG) argument of Ref. [l^. 
In this section, we further substantiate this argument. 
Starting from Eqs. ([2T]) and ([22|) . one can calculate 



the second-order correction to the bare vertex H^^^ = 
En Ke'""^ given by 



n(^)(P.,^) = ^in^EF„V^ 



(36) 



where we explicitly extracted the dependence on the up- 
per cutoff W . From Eq. (|36|) . we can immediately derive 
the following RG equations for the scale-dependent cou- 
plings rn(A = vfQY 



din A 2^ - 



(37) 



as in Ref. [23. This leads to the standard Cooper chan- 
nel renormalization. A direct derivation of these scaling 
equations can be found in Ref. [32|. At this lowest order, 
we obtain an infinite number of independent flow equa- 
tions, one for each angular momentum n. The integra- 
tion of these scaling equations directly leads to Eq. ([8]). 
These flow equations tell us that the couplings are 
marginally relevant in the infrared limit when the bare 
Tn are negative and marginally irrelevant otherwise. No- 
tice that at zero temperature, the running flow parameter 
A is replaced by the momentum vfQ in the Cooper chan- 
nel. The idea of the RG is to replace in the perturbative 
calculations of a momentum-dependent quantity the bare 
couplings Tn by their renormalized values. By doing so, 
we directly resum an infinite class of (ladder) diagrams. 

Fet us apply this reasoning now to the susceptibility 
diagrams and note that the first nonzero contribution to 
the linear-in-Q behavior of Xs{Q) appears in the second 
order in F^. For the particular example of ^Xs? renor- 
malization procedure has to be carried out independently 
for the two interaction lines, as illustrated by the series 
of diagrams in Fig. [71 For a given order of the inter- 
action ladder in the bottom (top) part of the diagram, 
one can perform the Cooper channel resummation of the 
top (bottom) interaction ladders to infinite order, as de- 
scribed in Sec. IIVBI or by using the RG equations. The 
fact that renormalized amplitudes F^ appear in the final 
results for the remaining diagrams (^Xi,2,4 is also clear 
from the RG argument, after insertion of particle-particle 
ladders as in Figs. [4]-[6l 

Finally, we note that the same series of diagrams that 
renormalizes the nonanalytic second-order contributions 
5xi 2 4 cilso contributes to the renormalization of the first- 
order diagrams displayed in Fig. [8] (notice that the first 
one is actually vanishing because of charge neutrality). 
As it is clear from the explicit calculation Sec. HVl the 
highest logarithmic powers, i.e., ex ihiVFQ /WY~^ at or- 
der n, renormalize Vm to F^ in the final expressions for 
Fig. [HI These first-order diagrams have an analytic de- 
pendence, at most Q^. Therefore, in agreement with the 
discussion in Sec. IIV A| the largest powers of the loga- 
rithms are not important for the linear dependence in 
Q and, in fact, they were already neglected to second 
order J 
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FIG. 8: First-order diagrams contributing to the spin suscep- 
tibility. These are renormahzed by the leading logarithmic 
terms of the higher order diagrams (see Figs. |4H6]). However, 
they do not produce a nonanalytic correction and can be ne- 
glected in the limit of small Q. 



VI. CONCLUSIONS 

In this paper we discussed the renormalization effects 
in the Cooper channel on the momentum-dependent spin 
susceptibility. The main result of the paper is given by 
Eq. (|35|) and shows that each harmonics gets renormal- 
ized independently. The derivation of the higher order 
corrections to the spin susceptibility was based on the 
second-order result, which we revisited through an inde- 
pendent direct calculation in the particle-particle chan- 
nel. Taking the angular dependence of the scattering po- 
tential explicitly into account, we verified that the main 
contribution indeed enters through forward- and back- 
scattering processes. At higher order, we found a simple 
and efficient way of resumming all the diagrams which 
contribute to the Cooper renormalization. We identi- 
fied the leading nonvanishing logarithm in each ladder 
and used this result in the second-order correction. This 
method saves a lot of effort and, in fact, makes the cal- 
culation possible. 

It was argued elsewhere that these renormalization ef- 
fects might underpin the nonmonotonic behavior of the 
electron spin susceptibility if the higher negative harmon- 
ics override the initially leading positive Fourier compo- 
nents. This would results in the negative slope of the spin 
susceptibility at small momenta or temperatures 



Other effects neglected here, as subleading logarithmic 
terms and nonperturbative contributions beyond the 
Cooper channel renormalization,^^ become relevant far 
away from the Kohn-Luttinger instability condition, but 
a systematic treatment in this regime is outside the scope 
of this work. Our results could be also extended to in- 
clude material-related issues such as disorder and spin- 
orbit coupling, which are possibly relevant in actual sam- 
ples. 

We also notice that final expression (|35|) parallels the 
temperature dependence discussed in Ref. |2l|, suggest- 
ing that the temperature and momentum dependence are 
qualitatively similar in two dimensions. This was already 
observed from the second-order calculation, in which a 
linear dependence both in Q and T is obtained. In our 
work we find that this correspondence continues to hold 
in the nonperturbative regime if the Cooper channel con- 
tributions are included. This conclusion is nontrivial and, 
in fact, does not hold for the three-dimensional case. 

The last remark, together with the experimental ob- 
servation of Ref. 24, supports the recent prediction that 
the ferromagnetic ordering of nuclear spins embedded 
in the two-dimensional electron gas is possible^i^ The 
ferromagnetic phase would be stabilized by the long- 
range Ruderman-Kittel-Kasuya-Yosida (RKKY) interac- 
tion, as determined by the nonanalytic corrections dis- 
cussed here. 
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APPENDIX A: DERIVATION OF THE LADDER DIAGRAMS 



To calculate the ladder diagram given by Eq. (|9|) we begin from the simultaneous change of all qi variables, 
Qi,n —Qi,iJ, —Pfi^ and expand the scattering potential into its Fourier components given by Eq. (|6]). Thus, 



mi...mn i=l 

(Al) 

where = ^ p^. We first evaluate the factors IiM{P^i) = ~ J -^^G{—q^)G{q^ + P^)e*^^'? appearing in the 
above formula. To this end, we integrate over the frequency Vtq and linearize the spectrum around the Fermi surface, 
^q+p ^ + vfP cos Oq. This rcquircs that 6q, and all the angles in Eq. ()Aip , are defined from the direction of P. 
We also use ^ = as integration variable, which gives 

n„(P„) = ,w,e'-'. hi eK + .,Pco.o )-e(-o ^ ^^^^ 



(27r)2 7o J 2^ ^vfP cos Oq-iftp 
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The unit step functions, + vpFcosOq) and 6(— determine the integration range in ^, which is —vpP cosOq to 
W and —W to 0, respectively, with W being the high energy cutoff. The energy integration yields 



m 



(27r) 



dOq e 



iTTSgnl^p + In 

2\ —vpF cos Uq — tilp 



In- 



2W 



vpP cos Oq — iftp 



(A3) 



with the sign term coming from the lower limit of the second integration, ln(— 21^ — iftp) — ln2W — ZTrsgnl^p. 

If we change variables in the integral over the first logarithm, 9q ^ 6q^ tt, we make it identical to the second one, 
except for the multiplicative term (—1)^ originating from g*^(^<?+^). Therefore, we find that 



nM(PM) 



(2 



dOq e 



iMOa 



sgmZp + m — h m ( 1 + i— — cos^ 



W 



(A4) 



where M is even and the factor of 2 in front of W has been absorbed into the cutoff. Writing the second logarithm 
as a series, ln(l + x) = — l)"^^;-, we can easily integrate term by term. 

The M = contribution is 



no(Pp) 



(27r) 



d0„ 



In- 



\np\ 
w 



^ n ( 

n—l 



-ivpP 

np 



cos 6. 



m 
2^ 



In- 



^ nUK 2np J 



w 



p 



(A5) 



where n is even in the primed sum. The summation of the series gives Eq. (pT]) shown in the main text 
For M 7^ (M even) we get 

' ' -ivpP 



m 



dOqC 



■ cos Oq J = 



m 
"2^ 



2>|M| 



(A6) 



Upon summation over even values of n, Eq. (p!2|) in the main text is obtained. We now consider Eq. ()Aip . and to 
simplify the notation we introduce summation indexes Mn-i = rrii — mi+i (for z = l,2,...n — 1) and m = —rrin- We 
also introduce m' = mi — = M^, which is even. We express the angles from the direction of p, i.e.. Op ^ —Op 
and Op' ^ — Op. Finally, we obtain 



m' ,Mi...Mn 



im Up — imb 



nMi(P^)...nM„_,(P^), 



(A7) 



which can be written as in Eq. (p!Q|) using definition (p!3|) . The special case of 11^^^ is given by Eq. ([2T]) and can be 
directly used for the derivation of RG equations in Sec. [Vl 



APPENDIX B: GREEN'S FUNCTIONS INTEGRATION FOR ^xi 



(n) 



We consider here the integration of the Green's functions appearing explicitly in Eqs. (p!8|) and ([23|) : 



hm 



1 



(27r)2 ^k/^^^k^ 27r (c^ + i^k)(cc; + ^^kO(^ + ^^k+Q)(^ - - ^^k-p) * 
The integration over co can be performed with the method of residues (we choose the lower half-plane contour); 



(Bl) 



A = im 



lim 



(27r)2 



e(a) 



0(Ck') 



(Ck' - ?k)(?k+Q - ik){^p + i^k + i^k-p) (■fk - Ck')(Ck+Q - Ck')(^P + iik' + «?k-p) 



e(a+Q) 



e(-a-p) 



(Ck - Ck+Q)('fk' - ^k+Q){^P + i^k+Q + i^k-p) i^P + iS,k + iS,k-p){^P + *Ck' + i^k-p){^P + iCk+Q + «Ck-p) 

(B2) 



Since the sum of all residues of the integrand in Eq. ()Bip in the entire complex plane is zero, we can subtract the 
same quantity ©(^k) from each numerator in Eq. (|B2p without affecting the result. This cancels the first term and, 
using lim^j^,^^j^ = ^(^k), we obtain 



A = Ai^A2^A3 =im 



f dOkO 
im / — — 
J (2^ 



dOkd^ 
? 



(^k - ^k+Q)(^P + ^^k + ^^k-p) 



e(a+Q) - ©(a) 



e(-a-p) - ©(a) 



(^k - ^k+Q)^(^P + ^^'k+Q + ^^k-p) (^P + ^^k + ^^k-p)^(^P + ^^k+Q + ^^k-p) 



(B3) 
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We can now perform the integration in d^i^. To this end we hnearize the energy spectrum near the Fermi surface, 
^k+Q ^ + ^fQ cos Ok and ^k-p ^ ^k — '^f^ cos Op^ which is a good approximation since Q^P <^ kp. We also define 
C = Vtp — ivpP cosOp (and the complex conjugate ( = ftp + ivpP cosOp). The first two integrals are vanishing, 



m 



27r 



dOk 



(vpQ cos Ok 



0, 



and 



m 



27r 



dOk 



C + ivpQ COS Ok 
{vpQ cos Ok^ \C ~ ivpQ cos Ok 



In 



0, 



since the integrands are odd with respect to cos Ok. 

For the remaining term we make use of the indefinite integral. 



/ 



z^-Z2 + {z^+2iO\n{^ 



{2i^ + zi)^{2i^ + Z2) 2i{zi- Z2)'^{2i^ + zi) 

which tends to zero for ^ ^ ±oo, and hence only the integration hmits at and vpP cosOp contribute; 



A = 



m 
8^ 



n 



1 



In 



In 



c 



{vpQ cos Ok)'^ I \C ^ ivpQ cos Ok J VC + ^'^fQ cos^/c 
The last term vanishes and the final integration in dOk can be done using 



ICP^fQcos Ok 



dOk 



I 

Jo 



o? cos^ 



dO 



27r 



1^1 - sgn(Re2;) j. 



which yields 
A = 



^ V(^p + ivpPcosOp)^ + {vpQ)^ ^ ^ 
2 - : sgnQp - : sgnl^p 



ftp + ivpP cosOp 



Qp — ivpP cosOp 



(B4) 



(B5) 



(B6) 



Notice that A is even in l^p. Furthermore, the last two terms give the same contribution upon integration J^^ dOp, 
and therefore Eq. ([19]) is obtained. 



APPENDIX C: CALCULATION OF ^Xi^^ 

In this appendix we consider the explicit evaluation of Eq. (|2Q|) . It receives contributions from all possible values 
of M appearing in 11*^^^ [see Eq. ([2T]) ]. We consider first M = for which Eq. ([22]) is useful. With the change of 
variables r = R{sm (j) -\- i cos (j) cos Op) we obtain 

(2) ''^ V ^^2 / 7/ / / jn I J 2/i V ' ^ J- \ V sm 0+2 cos 0COS 0P / \ W ) fr^^\ 

^Xi,o = -Tr^ — ^0 / d(^cos(^ dOp drr [I \ —3 , (CI) 

2ti^vf Jo Jo Jo \ r J (sm (p ^ i cos (p cos p^ 

where r^ax = '^^^{^^'^^ + i cos (p cos Op). Note, that it is necessary to introduce the upper cutoff in the integral over 
r, which is formally divergent. However, this upper limit turns out to be irrelevant for the nonanalytic correction. 
We start from the first contribution in the above equation, where the r integration gives 



Q drr'(l - y^^) = f [rL. - (l+rl^.f^' + 1] « ^(sin<A + icos0cos^p) + f . (C2) 

The term proportional to W, as in Ref. 8, is the dominant contribution to the spin susceptibility. However, it does not 
depend on Q and therefore is uninteresting for us. The term proportional to Q, from the lower integration limit, is 
the desired nonanalytic correction to the spin susceptibility and does not depend on (j) and Op. Therefore, the angular 
integration can be performed using^ 

7r/2 p27r In f 1+sin \ 

7, J TO V sin 0+i COS 0COS / 'T /r^o\ 

d(pcos(j) / dOp-— — ■ — — = --. (C3) 

Jo [smcp -\- 1 cos (j) cos Op j-^ 2 

The same analysis can be applied to the second term of Eq. ()Cip . which contains \n{vFQr /W). The integration in r 
gives a constant from the lower limit and the remaining angular integrations yield zero, as discussed in the main text. 
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Hence, the final result is 

t 

UttHf 



We now aim to calculate terms with M 7^ 0. By making use of Eq. (p!2|) for IImC^, and substituting again 
r = i?(sin (j) -\- i cos (j) cos Op) we obtain 



2|M|7r^VF Jo Jo Jo ^ r J \ zcoscj) J 



(sin + z cos cos ^p)^ 
(C5) 

The integral over r can be performed as before. The integration over Op yields 

/ 1 ^ = ^(^ + 3|M| sin(/> + 3 sin^ <\> - 1) ( / ) , (C6) 

Jq (sm0 + zcos^cost/p)'^^ V zcos0 / 

which can be obtained by standard contour integration in the complex plane [z = e~^^^). Finally, 

^XiM= «u7l^ E^M-„K r^'#cos0(i^^)'''''(M2 + 3|M|sin</>+3sin20-l) = -^^FM-nK, 

' Q\M\it^vf Jo ^ COS0 / Iztt^vf 

(C7) 

since the last integration gives a factor of |M|/2. Thus, the total second order correction is 

^xf^(Q) = ^^E'E2^M-„F„. (C8) 

M n 

Rewriting the double sum as ^ KnK, + Xlm n(~-'-)^^^^^^^' recover the two contributions as in Eq. ([T4j). 

APPENDIX D: SMALL Q LIMIT OF H^^^ 
We expand n*^^^ [see Eq. ([10])] in powers of Ilo{R^ (j)) as follows 

^(»)(i?,0,0p,^)«Fo"o...o^^l+ E'(^^Mo...o + ---^^oa^^^^ ■ (Di) 

The expression of no(i^, (^) is given by Eq. ([22]) . Therefore, the leading contribution in the above equation at small Q 
is from the first term since Hq"^ ^ In ^^^)^~^. However, the leading order does not contribute to the nonanalytic 
correction. Neglecting such constant terms, we can write the relevant subleading contribution in the following way: 

[no(i?, <t>)T-' = (n - 1) In ^)""' ^ In R{1 + sin </.) + ... = (n - 1) In !^)""'no(i?, </.) + .... (D2) 

Furthermore, by using Hq~^ ^ In ^^^)^~^, we can simplify Eq. (|Dip to the following form: 

n(")(i?, </., ^p, 0) = In !^)""'^' (y^, , + + . . . ^) Hm + . . . . (D3) 

M 

By making use of Eq. (p!3|) we have 

n-l 

^^Mo...o + v^,"m...o + ■ ■ • Ko...M = J2J2 yrvLM^'^''-''' m 

k j=l 

and therefore Eq. (jPSp is written explicitly as 



(D5) 
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We can now sum previous expression (|D5p over the index n > 2, which removes the constraint j -\- j' = n — 2. Hence, 
the last double summation factorizes in the product of two geometric series, that can be evaluated explicitly, and we 
obtain Eq. 
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